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Abstract 

Effect of considering stochasticity of Grassmannian coordinates in iV = 1 super- 
space was analyzed previously in the so called stochastic superspace framework of 
Kobakhidze et. al.. This naturally leads to the expected soft breaking terms of MSSM. 
The soft parameters like that of bilinear Higgs mixing, trilinear coupling as well as 
gaugino masses are proportional to a single mass parameter £, a measure of super- 
symmetry breaking arising out of the stochasticity. The soft parameters were given at 
a scale A that is in between the gauge coupling unification scale Mq and the Planck 
scale Mp. A non- vanishing trilinear coupling parameter is a natural outcome of the 
stochastic superspace framework. There is however no scalar field soft term at this 
scale. Confronted with tachyonic sleptons, scalar masses at the electroweak scale could 
only be radiatively generated if A > Mq. The model as such can not accommodate 
the lighter Higgs boson at 125 GeV. We extend the scope of this model in relation to 
the recent Higgs discovery, LHC limits of strongly interacting sparticle masses, WMAP 
data for dark matter, flavour physics constraints and XENON100 data by a simple and 
phenomenologically motivated addition of a non-vanishing scalar mass parameter at 
the high scale. In contrast to the previous analyses we consider A = Mq, thus avoiding 
any ambiguities of a post-GUT physics. 
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1 Introduction 



Low energy supersymmetry (SUSY) [1-4] has been one of the most promising candidate for 
a theory of fundamental particles and interactions, going beyond the Standard Model (SM); 
the so called BSM physics. The minimal extension of the SM including SUSY, namely the 
Minimal Supersymmetric Standard Model (MSSM), extends the particle spectrum of the SM 
by one additional Higgs doublet and the supersymmetric partners of all the SM particles - the 
sparticles. The Supersymmetric extension of the SM, provides a particularly elegant solution 
to the problem of stabilising the electroweak (EW) symmetry breaking scale against large 
radiative correction and keeps the Higgs 'naturally' light. In fact a very robust upper limit 
on the mass of the lightest Higgs boson is perhaps one of THE important prediction of this 
theory. Further, this upper limit is linked in an essential way to the values some of the SUSY 
breaking parameters in the theory. In addition in i?-parity conserving SUSY the lightest 
Supersymmetric Particle (LSP) emerges as the natural candidate for the Dark Matter (DM), 
existence of which has been proved beyond any doubt in the Astrophysical experiments. 
Search for evidence for realisation of this symmetry in nature, in the context of high energy 
collider experiments, precision measurements at the high intensity B-factories and in the 
DM detection experiments has therefore received enormous attention of particle physicists, 
perhaps only next to the Higgs boson. The recent observation of a boson with mass around 
~ 125 GeV at the Large Hadron Collider (LHC) [5] and the rather strong lower limits on the 
masses of sparticles which possess strong interactions that the LHC searches have yielded [6], 
necessitates careful studies of the MSSM in the context of all the recent low energy data. 
In these studies it is also very important to seek suitable guiding principles which could 
possibly reduce the associated large number of SUSY breaking parameters of MSSM. Thus 
looking for modes of specific SUSY breaking mechanisms that involve only a few input 
quantities given at a relevant scale can be useful. Here, the soft SUSY breaking parameters 
at the electroweak scale are found via renormalization group (RG) analyses. Apart from the 
simplicity of having a few parameters as input, such schemes create challenging balancing 
acts. On the one hand various soft breaking masses and couplings become correlated with 
one another in such schemes. On the other hand, overall one has to accommodate a large 
number of very stringent low energy constraints in a comprehensive model consisting of 
only a few parameters. A simple and well motivated example of a SUSY model is the 
minimal Supergravity (mSUGRA) [7]. Here SUSY is broken spontaneously in a hidden 
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sector and the breaking is communicated to the observable sector where MSSM resides 
via Planck mass suppressed Supergravity (SUGRA) interactions. The model involves soft- 
SUSY breaking parameters like i) universal gaugino mass parameter mi, ii) the universal 
scalar mass parameter m , hi) the universal trilinear coupling A , iv) the universal bilinear 
coupling B , all given at the gauge coupling unification scale. In addition to it one has 
the superpotential related higgsino mixing parameter /io with its associated sign parameter. 
The two radiative electroweak symmetry breaking (REWSB) conditions may then be used 
so as to replace B Q and /i with the Z-boson mass Mz and tan (3, the ratio of Higgs vacuum 
expectation values. Similar to mSUGRA one has other SUSY breaking scenarios like the 
gauge mediated SUSY breaking (GMSB) and models with anomaly mediated SUSY breaking 
(AMSB) etc [2,3]. Apart from direct collider physics data one has to satisfy constraints from 
flavor changing neutral current (FCNC) as well as flavour conserving phenomena like the 
anomalous magnetic moment of muon, constraints like electric dipole moments associated 
with CP violations or to check whether there is a proper amount of dark matter content in an 
R-parity conserving scenarios or proper neutrino masses in R-parity violating scenarios [1-3] . 
A single model is yet to be found that can adequately explain various stringent experimental 
results and at the same time possesses a sufficient degree of predictiveness. It is always 
however important to continue the quest of a simple model and check the degree of agreement 
with low energy constraints. 

In this work, we pursue a predictive theory of SUSY breaking by considering a field theory 
on a superspace where the Grassmannian coordinates are essentially fluctuating/stochastic 
[8-10]. We note that in a given SUSY breaking scenario, our limitation of knowing the 
actual mechanism of breaking SUSY is manifested in the soft parameters. Here, in stochastic 
superspace framework we assume that a manifestation of an unknown but a fundamental 
mechanism of SUSY breaking may effectively lead to stochasticity in the Grassmannian 
parameters of the superspace. With a suitably chosen probability distribution, this causes 
a given Kahler potential and a superpotential to lead to soft breaking terms that carry 
signatures of the stochasticity. As we will see, the SUSY breaking is parametrized by £ which 
is nothing but 1 /< 99 > , where the symbol < > refers to averaging over the Grassmannian 
coordinates. The other scale that is involved is A. Values of various soft parameters at this 
scale are the input parameters of the scheme. The values of the same at the electroweak 
scale are then obtained from from these input values, by using the the renormalization 
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group evolution. Considering the superpotential of MSSM, the soft terms obtained are 
readily recognized as the ones supplied by the externally added soft SUSY breaking terms 
of Constrained MSSM (CMSSM) [3], except that the model is unable to produce a scalar 
mass soft term [8]. Ref. [8] used A and £ as free parameters while analyzing the low energy 
signatures within MSSM. A was chosen between Mq to Mp, the scale of gauge coupling 
unification and the Planck mass scale respectively. The model as given in Ref. [8] and 
called as Stochastic Supersymmetric Model (SSM) is characterized by universal gaugino 
mass parameter mi, universal trilinear soft SUSY breaking parameter Aq, and universal 
bilinear soft SUSY breaking parameter Bq, all being related to £, the parameter related to 
SUSY breaking. We note that with the bilinear soft SUSY parameter being given, tan j3 
becomes a derived quantity. 

However, as already mentioned, in spite of the fact that the SSM generates soft SUSY 
breaking terms, it produces no scalar mass soft term. Scalar masses start from zero at the 
scale A and renormalization group evolution is used to generate scalar masses at the elec- 
troweak scale Mz- Scalar masses at Mz severely constrains the model because typically 
scalars are very light. In particular quite often sleptons turn to be the lightest supersym- 
metric particles (LSP) or even become tachyonic over a large part of the parameter space. 
This is partially ameliorated when one takes the high scale A to be larger than the gauge 
coupling unification scale Mq ~ 2 x 10 16 GeV. However, in spite of obtaining valid parameter 
space that would provide us with a lightest neutralino as a possible dark matter candidate 
in R-parity preserving framework, we must note that the low values that one obtains for 
the masses of the first two generation of squarks are hardly something of advantage in view 
of the constraints coming from FCNC as well as those from the LHC data [6]. On the 
other hand, SSM has a natural advantage of being associated with a non- vanishing trilinear 
coupling that is favourable to produce a relatively light spectra for a given value of Higgs 
boson mass m^. The recent announcement from CMS and ATLAS collaborations of the 
LHC experiment about the discovery of a Higgs like boson at ~ 125 GeV [5] thus makes 
this model potentially attractive. However, as explored in Ref. [10] SSM as such is unable to 
accommodate such a large rrih in spite of having a built-in feature of having a non-vanishing 
Aq. It could at most reach 116 GeV for rrih [9,10] and the constraint due to Br(B s — > 
as used in Ref. [10] was much less stringent in comparison to the same of present day [11]. 
Furthermore, it is also important to investigate the effect of the direct detection rate of dark 
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matter as constrained by the recent XENON100 data [12]. 

In this analysis we would like to give all the input parameters at the grand unification scale 
Mq, the scale at which Standard Model gauge group namely SU(3) x SU(2) x U(l) comes 
into existence. Any evolution above Mq would obviously demand choosing a suitable gauge 
group, a question that is not going to be addressed in this work. In this way we would like 
to avoid unknown issues arising out of a post-GUT [13] physics. However, we would rather 
try to meet the phenomenological demand of confronting the issue of sleptons becoming 
tachyonic or avoiding scalar masses to become light in general, in a minimal modification 
by considering an externally given scalar mass soft parameter mo as a manifestation of an 
additional origin of SUSY breaking. This would be useful to have adequately heavy first 
two generation of scalar masses so as to overcome the FCNC related constraints and LHC 
data [6] on squark masses. Additionally this will also be consistent with having the lighter 
Higgs boson mass (m^) to be in the vicinity of 125 GeV. We will henceforth denote the 
model as Mod-SSM. 

We may note that traditionally minimal versions of models of SUSY breaking have been 
extended for phenomenological reasons. It is also true that extending a minimal model often 
lowers predictiveness and may even cause partial dilution of the main motivations associated 
with the building of the model. For example, considering non-universal gaugino or scalar 
mass scenarios may be more suitable than CMSSM or mSUGRA so as to obtain a relatively 
lighter spectra in the context feasibility of exploring via LHC. Another example may be 
given in the context of the minimal AMSB model. As we know, a pure AMSB scenario [14] 
is associated with form invariance of the renormalization group equations (RGE) of scalar 
masses and absence of flavour violation. However, it produces tachyonic sleptons. In the 
minimal AMSB model [3, 15] one introduces an additional common mass parameter mo for 
all the scalars of the theory. This ameliorates the tachyonic slepton problem but it is true 
that we sacrifice the much cherished feature of form invariance and accept some degree of 
flavour violations at the end. We would like to explore a Non-minimal scenario of Stochastic 
Supersymmetry model, namely Mod-SSM in this spirit. We particularly keep in mind that 
the Stochastic Supersymmetry formalism may be used not only within MSSM framework 
but may be extended to superpotentials beyond that of MSSM [9]. The fact that the model 
with a minimal modification can easily accommodate the recent Higgs boson mass range by 
its generic feature of having a non- vanishing tri-linear coupling parameter makes it further 
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attractive. We believe that our approach of considering an additional SUSY breaking scalar 
mass term is justified for phenomenological reasons. 

Thus in Mod-SSM we consider the input of soft term parameters mi, B and A Q (all are 
either proportional to |£| or £* ) and a universal scalar mass soft parameter m , all being 
given at a suitable scale which we simply choose as the gauge coupling unification scale 
Mq considering the LEP data on gauge couplings could be a hint of the existence of Grand 
unification. As in Ref. [8] we would also restrict £ to be real, either positive or negative. 

2 Stochastic Grassmannian coordinates and SUSY break- 
ing 

As seen in Ref. [8] we consider an N = 1 superspace where the Grassmannian coordinates 9 
and 9 are taken to be stochastic in nature. One starts with identifying the terms involving 
superfields in the superpotential and the kinetic energy terms that could be used to construct 
the SUSY invariant Lagrangian density for a given model. Each term is then multiplied with 
a probability distribution function V(9, 9) and integrated over the Grassmannian coordinates 
appropriately. V(9, 9) can be expanded into terms involving 9 and 9 which obviously has a 
finite number of terms because of the Grassmannian nature of 9 and 9. One then imposes the 
normalisation condition / dWV(9,9) = 1 and vanishing of Lorentz non-scalar moments 
like < 9 >,< 9 >,< 99 >,< 9 2 9 > and < 99 2 >. The stochasticity parameter £ is defined 
as < 99 >= l/£*. Here, £ is a complex parameter with mass dimension unity. As computed 
in Ref. [8], and as worked out in this analysis explicitly in the Appendix, the above leads to 
the following Hermitian probability distribution. 

V(9, 9)\i? = H9, 0) = 1 + C{99) + i{99) + \i\\99){99). (1) 

For the simple case of a Wess-Zumino type of scenario [3] where the kinetic term is obtained 

1 1 

from and the superpotential being given as W = -m<3> 2 H — /i$ 3 , where $ is a chiral 

2 3 

superfield, one finds that the effect of stochasticity as described above leads to the following 
SUSY breaking term 

-L soft = l -Cm<p 2 + ^*h4> 3 + h.c. (2) 

Applying the stochasticity idea to the superpotential of MSSM, along with considering the 
effect on the gauge kinetic energy function, the above formalism leads to the following tree 
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level soft SUSY breaking parameters to be given at the high scale A: 

i) universal gaugino mass parameter mi = 

ii) universal trilinear soft parameter A = 2£*, 

iii) universal bilinear Higgs soft parameter Bq = £*. 

Recall that there is no scalar mass soft SUSY breaking term in SSM. 

For convenience we take £ to be a real positive number with an additional input sign(£). If 
we count the universal gaugino mass parameter mi as the independent parameter we have, 

A = sign(£)Ami,B = sign(£).2mi. (3) 

As has already been discussed before, we introduce a non-vanishing scalar mass parameter 
m and fix A at M^ 4 . Thus with the above extension, the input quantities for Stochastic 
SUSY model are: 

mi, m , sign(/i) and sign(£). 

We note that the model quite naturally is associated with non-vanishing trilinear soft break- 
ing terms. As we will see, this is quite interesting in view of the recent LHC announcement 
for the Higgs mass range centering around 125 GeV [5]. In this analysis, we will discuss only 
the case of £ < 0, because the other sign of £ does not produce a spectra compatible with 
the dark matter relic density constraint. 

The requirement of the REWSB then results in the following relations at the electroweak 
scale. 

2 _ 1 2 m\ D - m\ v tan 2 g E!-E 2 tan 2 /3 



and, 

sin 2(3 = 2Bfi/(m 2 HD + m 2 Hv + 2/i 2 + Ei + E 2 ) , (5) 

where Ej denote the one-loop corrections [16, 17]. Here, B refers to the value of bilinear 
Higgs coupling at the electroweak scale which has to be consistent with its given value Bq at 
Ma- Bq is determined via mi apart from a sign of the stochasticity parameter as mentioned 
before. Consequently, tan f3 is a derived quantity in the model. Bq at the scale Mq and B 



4 We note that a vanishing scalar mass parameter at a post-GUT scale, with RG evolution corresponding 
to an appropriate gauge group, would indeed generate non-vanishing scalar mass terms at the unification 
scale Mq [13]. 



7 



at the electroweak scale are connected via the following RGE written here at the one-loop 
level: 



where t = IniM^jQ 2 ^) with Q being the renormalization scale. Sj = a;/(47r) for i = 1,2,3 



running gaugino masses. Yj are the squared Yukawa couplings, e.g, Y t = y^/(47r) 2 where yt is 
the top Yukawa coupling. In this analysis the value of tan /3 is determined via Eqs.4, 5 and 
6 along with B = sign(^)2mi at the scale M G . We use SuSpect [18] for solving the RGEs 
and obtaining the spectra. The code takes tan /3 as an input quantity. Hence, we implement 
a self-consistent method of solution that starts from a guess value of tan /3 resulting into a 
B(Mg) that in general would not agree with the input of B . Use of a Newton-Raphson root 
finding scheme ensures a fast convergence toward the correct value of tan/3 when B(Mq) 
matches with the input of Bq. Here we stress that we do not encounter any parameter point 
with multiple values of tan /3 in our analysis 5 . 



The fact that the model has tan (3 as a derived quantity necessitates studying the behaviour 
of the evolution of the bilinear Higgs parameter B. Fig.l shows the evolution of a few rela- 
vant couplings for a specimen input of mi/2 = 600 GeV, mo = 2 TeV and fi > in Mod-SSM. 
For £ < 0, both Bq as well as Aq are negative namely Bq = — 2m 1 / 2 and Aq = — 4m!/ 2 . For 
a valid parameter point within the model with /i > we require B > 0, a necessity to have 
a positive sin 2/3 from Eq.5 6 . We note that the denominator in the right hand side of Eq.5 is 
the square of pseudoscalar Higgs mass which needs to be positive. The fact that B is origi- 
nally negative at Mq and has to change to a positive value at Mz puts a strong constraint 
on the parameter space of the model. Numerically, this results into tan/3 assuming large 
values. In regard to the evolution of A-parameters we defer our discussion until Fig. 3. 

Fig. 2 shows scatter plot of parameter points in tan /3 — mi plane that satisfy the REWSB 

5 See Ref. [19] for such a general possibility in REWSB where Bo is given as an input. 

6 tan/3 which is the ratio of two vacuum expectation values is positive. Hence sin 2/3 — n+tan^fl) ^ s a ^ so 



dB 




(6) 




3 Results 



positive. 
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Figure 1: Evolution of a few relavant couplings for a specimen input of m\/2 = 600 GeV, itlq = 
2 TeV and \x > 0. With £ < 0, one has Bq = —2m\j2 and Aq = — 4m 1 / 2 - For a valid parameter 
point within the model with /i > Owe require B = B{Mz) > 0, a necessity to have a positive sin 2/3 
from Eq.5. 



constraints of Eqs.4 and 5. As m is varied up to 7 TeV and mi up to 2 TeV, tan (3 is 
seen to have a range 32 to 48. The spread of tan j3 for a given mi arises from variation of 
m . This spread decreases as mi becomes larger because of increased sensitivity to satisfy 
REWSB. It is worth mentioning that the range of valid tan (3 is much larger for the case 
of £ > where B stays positive throughout the range from Mq to Mz- This is unlike the 
case of £ < under discussion where B is negative at Mq and necessarily has to become 
positive at the electroweak scale, thus adding stringency to tan /3 in its range. However, as 
already mentioned, we will not discuss the case of £ > further because of the resulting 
overabundance of dark matter for the entire parameter space for this sign of £. 

We will study now the effect of low energy constraints particularly in the context of the 
recent discovery of Higgs-like boson [5]. Fig.3 shows the result in the mi — m plane for 
£ < 0. A sufficiently non-vanishing A t (see for example Fig.l) helps in producing a large 



9 



60 
55 
50 
45 
40 
35 

CO. 

§ 30 

-4—* 

25 
20 
15 
10 
5 



i;<o 




_L 



_L 



_L 



_L 



200 400 600 800 1000 1200 1400 1600 

m 1/2 (GeV) 



1800 2000 



Figure 2: Scatter plot of parameter points in mi — tan /3 plane when mo and mi are scanned up to 

2 2 

7 TeV and 2 TeV respectively. Here we only consider the validity of REWSB constraints of Eqs.4 
and 5. The values of tan/3 that satisfies the REWSB constraint varies from 32 to 48. The spread 
of tan/3 for a given mi arises from variation of mo- This spread decreases as mi becomes larger 

2 2 

because of increased sensitivity to satisfy REWSB. 



loop correction to the lighter CP-even Higgs boson h. This can be understood by looking 
at the expression for the dominant part of loop correction to the Higgs boson mass coming 
from the top-stop sector [20-22] . 



Am 2 h 



3m t 



^ 4 



2ir 2 v 2 sin f3 



log 



Ml 
ra 2 



2M| 



1 - 



6M| yj 



(7) 



Here Ms = yffn^ rnf 2 , X t = A t — /xcot 0, v = 246 GeV and m t is the running top quark 
mass that also takes into account QCD and electroweak corrections. The loop correction is 
maximized if X t = y/6Ms- Clearly, a non- vanishing Aq can be useful to increase Am| so that 
rrih reaches the LHC specified zone without a need to push up the average sparticle mass scale 
by a large amount. As particularly mentioned in Ref [8], the SSM is additionally attractive in 
this context since it naturally possesses a non- vanishing and large \A \. We further note that 
in this model the lighter Higgs boson h has couplings similar to those in the Standard Model, 
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because the CP-odd Higgs boson mass (tha) is in the decoupling zone [23]. The ATLAS and 
CMS results for the possible Higgs boson masses are 126.0 ± 0.4 (stat) ± 0.4(syst) GeV and 
125.3 ± 0.4 (stat) ± 0.5(syst) GeV respectively [5]. In regard to MSSM light Higgs boson 
mass we note that there is about a 3 GeV uncertainty arising out of uncertainties in the top 
quark mass, renormalization scheme, as well as scale dependence and uncertainties in higher 
order loop corrections up to 3- loop [24-28] . Hence, in this analysis we consider the following 
limits for m^. Two blue lines show the corresponding limits in Fig.3. 

123 GeV < m h < 127 GeV. (8) 

We will outline the other relevant limits used in this analysis. A SUSY model parameter 
space finds a strong constraint from Br(b — > 57). In SM, the principal contribution that 
almost saturates the experimental value comes from the loop comprising of top-quark and 
W-boson [29]. In MSSM, principal contributions arise from loops containing top-quark 
and charged Higgs bosons, and the same containing top squarks and charginos [30]. The 
chargino loop contributions are proportional to A t [i and this may cause cancellations or 
enhancements between the principal terms of the MSSM contribution depending on the 
sign of A t [i. Similar to mSUGRA, both SSM and Mod-SSM with £ < also typically have 
At < 0. With > 0, this primarily means cancellation between the chargino and the charged 
Higgs contributions. This leads to a valid region for Br(b — > S7) for a larger sparticle mass 
scale compared to the case of /1 < 0. Br{b — > sj) constraint thus favours the positive 
sign of /i by allowing larger areas of parameter space. We consider the experimental value 
Br(b — > 57) = (355 ± 24 ± 9) x 10 -6 [31]. This results into the following 3a level zone as 
used in this analysis. 

2.78 x 10~ 4 < Br(b s 7 ) < 4.32 x 10" 4 . (9) 

The above constraint is displayed as a maroon dot-dashed line. The left region of this line 
would be a discarded zone. Next, the fact that the Stochastic Model with £ < selects 
appreciably large values for tan/3 necessitates checking the B s — > fi + n~ limit. This is 
required because B s — > increases with tan (3 as tan 6 (3 and decreases with increase in 

viiA-, the mass of pseudoscalar Higgs boson as m^ 4 [32]. We have used the recent experimental 
limit [11] at 95% confidence level as follows. 

Br(B s < 4.5 x 10~ 9 . (10) 
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Figure 3: Constraints shown in the mi — m plane for fj, > and £ < for Mod-SSM. A and B 

2 

in the model satisfy Aq = sign(^)Ami and Bq = sign(^).2mi . tan/3 becomes a derived quantity 

2 2 

that varies between 32 to 48. The Higgs boson limits are shown as two solid blue lines. Br{b — > S7) 
limit is shown as maroon dot-dashed line. The lower part corresponds to discarded region via Eq.9 
where the branching ratio goes below the lower limit of the constraint. Br(B s — > fi + fi~) limit is 
shown as a brown solid line of which the lower region exceeds the limit of Eq.10. The top green 
region (I) corresponds to discarded zone via REWSB. The bottom blue-green region (III) refers 
to the zone where stau becomes LSP or tachyonic. The gray region (II) has discontinuous patches 
of valid parameter zones, the details of which are mentioned in the text. Red points/areas falling 
in region-II satisfy the WMAP-7 data only for the upper limit of Eq.12. Typically, the red points 

bordering region-Ill and some part of the extreme left red points (for very small mi) satisfy both 

2 

the upper and the lower limits of Eq.12. 



Br(B s — y /i + /i~) limit is shown as a brown solid line going across Fig.3, of which the lower 
region exceeds the above limit. We also compute Br(B + — > t + v t ) in this analysis. The SUSY 
contribution to Br(B + — > t + u t ) is typically effective for large tan (3 and small charged Higgs 
boson mass scenarios [33]. We consider the experimental data from BABAR [34] which is 
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Br(B + -»■ r + v T ) = 1.83j£i(stat.) ± .24(syst.)) x 10" 4 . This has to be compared with 
the SM prediction (1.18 ± 0.16) x 10 -4 [34]. Adding the errors in the experimental result in 
quadrature we obtain the following : 



Br(B + — > r + u Tl r 

R(b^t+v t ) = D \ D+ t ; f SY = 1-55 ± 0.34. 11 

Br(B+ r+u T ) SM 



This translates into 0.53 < R(b+^t+u t ) < 2.57 at 3cr. Here Br(B + — > i~ + v T ) SUSY denotes the 
branching ratio in a SUSY framework that also includes the SM contribution. In general, we 
find that the model parameter space of Mod-SSM is not constrained by B + — > t + u t since 
charged Higgs bosons are sufficiently heavy. 

We have not however included the constraint from muon g — 2 in this analysis considering the 
tension arising out of large deviation from SM value, uncertainty in hadronic contribution 
evaluations and accommodating SUSY models in view of the LHC sparticle mass lower 
limits [35]. 

We now explore the cosmological constraint for neutralino dark matter relic density [36]. At 
3cr, the WMAP-7 data [37] is considered as shown below. 

0.094 < tt~oh 2 < 0.128. (12) 

The conclusions in regard to the relic density constraint is additionally found to be sensitive 
on the top-quark mass in this model. We divide the dark matter analysis into two parts 
depending on a) the top-quark pole mass set at 173.3 GeV and b) using a spread of top-quark 
pole mass within its range m t = 173.3 ± 2.8 GeV following the result of the recent analysis 
performed in Ref. [38]. In this context we note that the experimental value as measured in 
CDF and DO collaborations of Tevatron is mT v = 173.2 ± 0.9 GeV [39] 7 . 



3.1 Analysis with m t = 173.3 GeV: Underabundant LSP 

The lightest neutralino, the lightest stable particle (LSP) of the model is typically highly 
bino dominated except in a few regions where the higgsino mixing parameter fi turns out 
to be small. The parameter points within all the white region in Fig. 3 have bino-dominated 
LSP. At this point we note that the implementation of the REWSB conditions as manifest in 
Eqs.4 and 5 has to be done keeping in mind i) positivity of sin 2f3 and ii) positivity of fi 2 and 

7 For bottom quark mass we have used mt MS (m(,) = 4.19GeV. 
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and iii) the requirement of Bq, Aq being related to mi as given by Eq. 3, as well as requiring 
that lighter chargino mass lower limit is respected. All these requirements lead the top green 
shaded region (labeled as "I") to be a discarded zone. On the other hand, the gray shaded 
region (II) has discontinuous zones of valid parameter points, shown in red. The red points 
have considerably small values of \x thus giving the LSP a large degree of higgsino mixing. 
Apart from the red points, there are no solutions in the gray areas of region II. The need to 
satisfy the point (i) as explained above along with the requirement to satisfy the condition 
(iii) (implemented via a Newton- Raphson method of finding the correct tan/3) stringently 
negates existence of solution zones within region II. As a result, either there are solutions 
with appreciably small fi or no solution at all within this region. This on the other hand 
leads to a large amount of Xi ~~ Xi coannihilation. Such degrees of coannihilations indeed 
cause LSP to be only a subdominant component of DM. Thus, in this part of the analysis 
we consider the possibility of an underabundant dark matter candidate and ignore the lower 
limit of the WMAP-7 data. Typically we see that the relic density falls below the lower limit 
of Eq.12 in region II by an order of magnitude. Such underabundant LSP scenarios have 
been discussed in several works [40]. 

The lower shaded region III is disallowed as the mass of the stau (ti) turns negative or it is 
the lightest supersymmetric particle (LSP). Typically the red strip near region III refers to 
the LSP-stau coannihilation 8 zone where the relic density can be consistent with both the 
upper and lower limits of Eq.12. Quite naturally, the coannihilation may be stronger and 
this would additionally produce some underabundant DM points for this zone. Finally, the 
leftmost red region (with very small mi) satisfying WMAP-7 data is discarded by all other 
constraints. We note that only a small region satisfying the Higgs mass bound is discarded 
via Br(b — > sj). On the other hand, the Higgs mass bound line of 123 GeV supersedes 
the constraint imposed by the recent data on inclusive search for SUSY by the ATLAS 
experiment [6]. The most potent constraint to eliminate a large region of parameter space 
with m up to 2 TeV or so is due to Br(B s — > fi + fi~) data (Eq.10). The constraint is effective 
simply because of the large values of tan /3 involved in the model. The discarded parameter 
space via the above constraint includes the stau-coannihilation zone for dark matter. 

We will now describe the spin-independent direct detection scattering cross section for 

8 See for example works in Ref. [41] for various annihilation processes in relation to SUSY parameter space 
in general. 
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scattering of LSP with proton. The scalar cross section depends on t-channel Higgs exchange 
diagrams and s-channel squark diagrams. Unless the squark masses are close to that of LSP 
the Higgs exchange diagrams dominate [42]. We note that for the cases of parameter points 
with < (QcDMh 2 ) min , where (flcDMh 2 ) min refers to the lower limit of Eq.12, one must 
appropriately include the fraction of local DM density contributed by the specific candidate 
of DM under discussion while evaluating the event rate. 




200 400 600 800 1000 



m (GeV) 

Figure 4: Scaled spin-independent Xi~ P scattering cross section vs LSP mass. The scaling factor 
is given as C = ^^h 2 /(^lcDMh 2 ) min , where (QcDMh 2 ) min refers to the lower limit of Eq.12. 



This translates into multiplying <r|~o for such underabundant scenarios by p x /po- Here p x is 
the actual DM density contributed by the specific DM candidate contributing to po where the 
latter is the local dark matter density. We thus use p x = p C where ( = fl~h 2 / (QcDMh 2 ) min . 
On the other hand ( is simply 1 for abundant or over-abundant dark matter cases. Thus, we 
conveniently define ( = min{l, Q~h 2 / (QcDMh 2 ) min } [43]. Fig.4 shows the rescaled cross sec- 
tion as computed via micrOMEGAs version 2.4 [44]. We wish to emphasize that while some 
region of parameter space where LSP typically has a large degree of higgsino component is 
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eliminated via XENON100 data as announced in the summer of 2012 [12], a large section 
of parameter space remains to be explored via future direct detection of DM experiments. 
This consists of both types of coannihilation zones namely the chargino as well as the stau 
coannihilation zones. We must also keep in mind the issue of theoretical uncertainties, par- 
ticularly the hadronic uncertainties in evaluating 0"p~o- The strangeness content of nucleon 
finds a large reduction in the evaluation of relevant couplings via lattice calculations [45]. 
This is not incorporated in our computation while using micrOMEGAs to calculate the cross 
section. Thus, the above itself will cause a reduction of a s ~ Q by almost an order of magnitude. 
There is also an appreciable amount of uncertainty of the local dark matter density [46] . All 
these points need to be kept in mind while evaluating the implications of Fig. 4 for SUSY 
models. 



3.2 Analysis with m t = 173.3 ± 2.8 GeV: LSP of right abundance 

It is to be noted that in the gray area (shown as region II) of Fig. 3 the existence of a valid 
solution depends very critically on the parameters of the model. It was found that either i) 
we obtain a very small \x (barely satisfying the lighter chargino mass lower limit) in the gray 
area that would only provide us with extreme coannihilation between lighter chargino and 
LSP leading to underabundance of DM or ii) we find no valid solution at all. The sensitivity 
arises from stringency of satisfying REWSB on the parameter space for this region. In 
other words, for a given mi, a small change in assumed mo for a valid parameter point 
would produce a small change in tan /3 in the white region producing most probably another 
valid parameter point. However, the change may not be allowed via REWSB, particularly 
Eq.5 if the parameter point is considered in the gray region. Eq.5 means sin 2/3 needs to 
be a positive quantity less than unity. This typically becomes a severe constraint even if 
condition of satisfying the lighter chargino mass lower limit is met. With the top quark mass 
having a strong influence on REWSB it may be useful to investigate whether varying m t 
may extract newer valid points within the gray region that would have a suitable /i so as to 
satisfy a well-tempered [47] LSP situation. This will then have the right abundance of dark 
matter satisfying both the lower and the upper limit of DM of Eq.12. 

The present experimental data on top quark mass reads m e * v = 173.2 ± 0.9 GeV [39]. 
Recently Ref. [38] predicted the pole mass of top quark to be mf ole = 173.3 ± 2.8 GeV. 
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The analysis used the next-to-next-to leading (NNLO) in QCD prediction of the inclusive 
pp — >■ ti + X cross section and the Tevatron and LHC data of the the same cross section. 
The comparison between the experimental and the theoretical results helped extracting the 
top quark mass in the modified minimal subtraction (MS) scheme. This was then used 
to compute the pole mass m^ ole . We now extend our analysis by investigating the effect 
of varying top quark pole mass within the above range (mj? olc = 173.3 ± 2.8 GeV) on the 
solution space of the model. Indeed, we will see that even a variation of 0.9 GeV, the range 
of experimental error would be enough to have substantial effect on the conclusions. 
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Figure 5: Relic density satisfied points shown in red that fall within the lower and upper limits 
of WMAP-7 data for the neutralino relic density, when mt is varied by 2.8 GeV on either side of 
173.3 GeV. Here mi and mo are scanned up to 2 TeV and 7 TeV respectively for a > and £ < 0. 

2 

The two blue lines are the WMAP-7 limits of Eq.12. The value of mt = 173.3 GeV as used in Fig. 3 
visibly falls in a disfavoured zone in the context of obtaining the correct relic density, particularly 
toward the lower limit of of Eq.12. 



Fig. 5 shows the scattered points that satisfy the WMAP-7 given range of relic density for 
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the aforesaid variation of m t = m^ oc . Here, mi and m are varied up to 2 TeV and 7 TeV 
respectively. We see from this figure that the central value of m t (=173.3) GeV as considered 
in Fig. 3 is indeed the one which has the least amount of possibility to satisfy the WMAP-7 
data. The occurence of points which satisfies WMAP-7 data is particularly rare around this 
value of m t , near the lower part of the limit of relic density . We already found that the 
gray region (region-II) of Fig. 3 is a sensitive zone because of REWSB, where /i can be quite 
small. For such small values of \i one can only expect a large degree coannihilation 
which results into very small relic density. The latter goes below the lower limit of Eq.12. 
Hence no red point exists near the bottom blue line of Fig. 5 for this value of m t . On the 
contrary, a value of m t less than 1 GeV from the central value, which is only within the 
experimental error of m£ xp , would cause to have an LSP with correct abundance for DM. 
The most favoured zone for m t however would be from 171 to 172 GeV for satisfying the relic 
density limits. In fact, a reduced degree of sensitivity for satisfying REWSB is the reason 
of obtaining a well-tempered LSP while considering a top quark mass little away from the 
central value. 

Fig. 6 shows the effect of scanning the top quark mass on mi — m plane. Region I shown 

2 

in green is a disallowed area of parameter space via REWSB similar to Fig.3, except that it 
is now an invalid region for all values of m t within its limit. Thus this region is a smaller 
in extension than the corresponding region of Fig.3. The region II is disallowed because of 
stau turning to be least massive or being tachyonic. The constraints from Br{b — > S7) and 
Br(B s — > are as shown. The lines denote the boundary of purely discarded zones 

irrespective of variation of m t within its range. The blue line for = 124 GeV means that 
m/j < 124 GeV for all the region left of the line irrespective of the value of m t within the 
range. The ATLAS specified limit [6] for squarks also falls well within this left zone. The 
neutralino relic density satisfied areas are below region I and above region II. The red points 
satisfy both the limits of the WMAP-7 data, thus correspond to having the right degree of 
abundance of DM. On the other hand we have also shown blue-green points that only satisfy 
the upper limit of the WMAP-7 data. In this part of the analysis we consider LSP to have 
the correct abundance so as to be a unique candidate for DM. 

Finally we show the effect of varying m t on the spin-independent LSP-proton scattering 
cross-section in Fig. 7. Only the WMAP-7 satisfied points are shown (in red). Considering 
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Figure 6: Scattered points in mi — tuq plane for fi > and £ < when top quark mass is varied 

2 

by 2.8 GeV on either side of 173.3 GeV. Region I is same as that in Fig. 3 except that it is now a 
discarded zone for all values of mt within its limit. The region II is discarded via stau turning the 
LSP or turning itself tachyonic for all mt- The constraints from Br(b — > s-y) and Br(B s — > 
are as shown. The lines show the boundary of purely discarded zones irrespective of variation of 
mt within its range. There is no region toward the left of the blue line labeled as = 124 GeV for 
which rrih may become larger than 124 GeV irrespective of values of mt- Three benchmark points 
A,B and C are shown corresponding to Table 1. 

an order of magnitude of uncertainty (reduction) primarily because of the issue of strangeness 
content of nucleon as well as astrophysical uncertainties as mentioned before, we believe that 
the recent XENON100 data still can accommodate Mod-SSM even while considering the LSP 
as a unique candidate of DM. 

Table 1 shows three benchmark points of the model. The top quark mass m t is as shown 
for each of the cases. Points A and B correspond to values of m t which are entirely within 
the experimental error. The points A and C correspond to the upper region of Fig. 6 and 
these two points correspond to hyperbolic branch (HB)/focus point (FP) zone [48,49]. These 
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Figure 7: Spin-independent LSP-proton scattering cross-section vs LSP mass when mt is scanned 
along with mi and mo for \l > and £ < 0. Only WMAP-7 satisfied points (for both the lower 
and the upper limit) are shown along with the XENON100 exclusion limit. 



points are associated with a large degree of Xi ~ xt coannihilation. The degree of agreement 
between the desired value of Bo and the one obtained via Newton-Raphson iteration may be 
seen from the 5th and the 6th rows. We allowed a maximum deviation of 1 percent in the 
iterative procedure while scanning the parameter space. Majority of the points are consistent 
within tenth of a percent in this regard. The charginos and the neutralinos are relatively 
lighter for points A and C in comparison to those in point B. The scalar particles on the other 
hand are relatively lighter for point-B. Undoubtedly, like many SUSY models analyzed or 
reanalyzed after the Higgs boson discovery, the spectra is on the heavier side. The fact that 
Aq being non-vanishing and adequately large helps in reducing the average sparticle mass 
to a great extent compared to vanishing Aq scenarios satisfying the current lighter Higgs 
boson limit, rrih for point-C however goes below the assumed limit of Eq.8. We however still 
believe that it is within an acceptable zone considering the various uncertainties to compute 
rrih as mentioned before. Point-B has a slightly larger value of Br(B s — > than the 
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considered limit. Point A and C have larger spin-independent scattering cross section 
than the XENON100 limit. But, we believe this is within acceptable limit considering the 
existing uncertainties arising out of strangeness content of nucleon as well as those from 
astrophysical origins particularly from local DM density. Finally, it is also possible to satisfy 
all the limits in full subject to a five to ten percent heavier spectra and/or considering a 
multi-component DM scenario. 
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Parameter 


A 

A 


B 


C 


m t 


173.10 


173.87 


171.58 


m 1/2 


838.78 


1239.16 


579.69 


m 


6123.75 


1817.69 


4200.55 


(Aq = -4m 1/2 ) 


-3355.13 


-4956.64 


-2318.75 


(B = -2m 1/2 ) 


-1677.56 


-2478.32 


-1159.37 


i?o (as output) 


-1683.56 


-2478.32 


-1160.27 


tan /3 (as output) 


45.86 


40.92 


45.11 


sgn(n) 


1 


1 


1 


fi 


403.86 


2508.85 


310.43 


m-g 


2145.53 


2727.64 


1525.80 


niu L 


6247.84 


2994.87 


4292.70 




3758.76, 4376.60 


1333.10, 2078.60 


2587.56, 3026.20 


™~ bl ,m~ b2 


4397.10, 4886.58 


2054.58, 2339.25 


3037.22, 3381.67 




6119.53, 6119.05 


1983.72, 1982.21 


4197.61, 4196.89 


m fl , m y ~ T 


4750.43, 5482.91 


549.62, 1536.41 


3281.31, 3770.25 




406.25, 741.03 


1038.00, 2491.42 


304.84, 518.70 


X\ 1 A2 


356.62, 417.02 


548.95, 1038.00 


241.80, 316.70 


A3 ' A4 


424.40, 741.06 


2489.56, 2491.14 


322.35, 518.80 


m A ,m H ± 


2573.69, 2573.69 


1967.83, 1967.50 


1846.04, 1846.04 


m h 


124.42 


126.55 


121.57 




0.1105 


0.1002 


0.1106 


BF(b -»■ s 7 ) 


3.23 x 10" 4 


2.96 x 10" 4 


3.15 x KT 4 


BF(B S -> 


2.98 x 10^ 9 


5.23 x 10~ 9 


2.89 x 10~ 9 




0.98 


0.98 


0.97 




5.78 x 10~ n 


1.65 x 10" 10 


1.22 x 10~ 10 


a 57 


3.05 x 10~ 8 


1.04 x 1Q- 11 


2.64 x 10~ 8 



Table 1: Spectra of three specimen parameter points A, B and C as shown in Fig. 6. Results 

with marginal deviation from the assumed limits are shown in italics (red). Muon (7 — 2 is 

not imposed as a constraint in this analysis. -Bo = B(Mq) has two entries. The first one is 

the desired value while the second one is the value obtained with a suitable tan /3 as found 

from the Newton-Raphson root finding scheme. See text for further details. 
F 6 22 



4 Conclusion 



In Stochastic Superspace Model (SSM) [8] one assumes that a manifestation of an unknown 
but a fundamental mechanism of SUSY breaking may effectively lead to stochasticity in 
the Grassmannian parameters of the superspace. With a suitable probability distribution 
decided out of physical requirements, stochasticity in Grassmannian coordinates for a given 
Kahler potential and a superpotential may lead to well known soft breaking terms. When 
applied to the superpotential of MSSM, the model leads to soft breaking terms like the 
bilinear Higgs coupling term, a trilinear soft term as well as a gaugino mass term, all related 
to a parameter £, the scale of SUSY breaking. The other scale considered in the model of 
Ref. [8] at which the input quantities are given is A, where the latter can assume a value 
between the gauge coupling unification scale Mq and the Planck mass scale Mp. There is 
an absence of a scalar mass soft term at A in the original model that leads to stau turning 
to be the LSP or even turning itself tachyonic if A is chosen as Mq. This is only partially 
ameliorated when A is above Mq. However, the model because of its non- vanishing trilinear 
parameter is potentially accommodative to have a larger lighter Higgs boson mass via large 
stop scalar mixing. As recently shown in Ref. [10] the model in spite of its nice feature of 
bringing out the desired soft breaking terms of MSSM is not able to produce rrih above 116 
GeV and it has an LSP which is only a subdominant component of DM. 

In this work a minimal modification (referred as Mod-SSM) is made by allowing a non- 
vanishing single scalar mass parameter mo as an explicit soft-breaking term. Phenomeno- 
logically the above addition is similar to what was considered in minimal AMSB while 
confronting the issue of sleptons turning tachyonic in a pure AMSB framework. The mod- 
ified model successfully accommodates the lighter Higgs boson mass near 125 GeV. Addi- 
tionally, it can accommodate the stringent constraints from the dark matter relic density, 
Br(B s — > Br(b — > 57) and XENON100 data on direct detection of dark matter. A 

variation of top quark mass within its allowed range is included in the analysis and this shows 
the LSP to be a suitable candidate for dark matter satisfying both the limits of WMAP-7 
data. 
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5 Appendix 

Using Ref. [8] we consider the following Hermitian probability distribution : 

V(6, 9) = A + 6 a V a + d H d + e a 9 a B + 9~J & C + 9 a a\ p 9% 

+9 a 9j a A. & + ej"0 a ^a + 6 a 8jJ«D (13) 

Here A, B, C, D and are complex numbers. S, A, and S are Grassmann numbers. 
In order to arrive at the results of Ref. [8], we use the following [3]: d 2 9 = — \d9 a d9 a) 

d 2 e = -lde & de & , d 4 e = d 2 ed 2 e, jd 2 e (99) = 1, jd 2 e (99) = 1, jd 2 e = jd 2 e = o and 
/ d 2 e e a = j d 2 e 9 & = o. 

Normalisation: 

First, V(9,6) should satisfy the normalisation condition / d 2 9d 2 9 V(9,9) = 1. All the terms 
except the one with the coefficient D vanishes in J d 2 9d 2 9 V(9,9). Thus D — 1. 
Vanishing fermionic moments: 

Next, we require vanishing of moments of fermionic type because of the requirement of 
Lorentz m _ T he faot < * >-j Sm ^ S) - means E „ . 0. S^, 
< 9p >= means A a = 0, and < 6^% >= leads to = 0. Finally, < 9 2 9^ >= gives 

= and < 9 P 9 2 >= gives * a = 0. 
Bosonic moments: 
We now compute the bosonic moments. 
<99>= j d 2 9d 2 9 e p e p V{e,9) = J d 2 9d 2 9 (99)(99)C = C. 
Similarly, < 99 >= B. Calling B = l/£ one has B = C* = l/£. 

Furthermore,< 9999 >= A. The fact that < 9999 >=< 99 >< 99 > leads to A = 1/|£| 2 . 
Thus we find the following Hermitian probability measure for the stochastic Grassmann 
variables. 

V(9, 9)\Z\ 2 = V(9, 9) = 1 + C{09) + i{99) + \i\ 2 {99){99). (14) 
We consider Wess-Zumino model with a single chiral superfield $. $ has the following 
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expansion [3]. 



(x) - iBa^Bd^x) - -9 2 9 2 d tl d li (f)(x) + V29if>{x) 



+^e 2 d fl ^(x)a> i e + e 2 F(x) 

v2 



(15) 



Correspondingly for $t we have, 

$t = ( h*(x)+i9a fl 9d u (f)*(x) 



]e 2 9 2 d fl dy*(x) + V29^(x) 



— r 9 2 9a»d^(x) + 9 2 F*(x). 
V 2 



(16) 



The kinetic term of the Lagrangian £ is 



$t$ 



. One finds, 



D 



$t$ = |0| 2 + y2^0* + V29ip(p + 9 2 (I)*F + 9 2 F*^) + 29ip9ip 
+iV29 2 9a (1 ilj[d fl }(t)* + \p29 2 9^F - 2i9a^9(j)*[d^ 
+iV29 2 9a^[d„](f) + y/29 2 9i/;F* 



+9 2 9 2 ( F*F 



(17) 



Here, X[<9 M ]F = ~ (Xc^F — F<9^X). Upon using vanishing of appropriate surface terms, the 
D-term in particular reads: 

$t$] d = F*F + 8^*8^ + i - d^a^) . (18) 

Next, we consider the superpotential W, given by W = im$ 2 + |/i$ 3 . One has, 

$ 2 = (f) 2 + 2V29ip(j) + 9 2 (2(j)F -tpi/j), and 

$ 3 = (j? + 3V29ip(j) 2 + 39 2 (F(j) 2 (19) 

Thus, 

W = Qm0 2 + \h(jy^j + V294)(m4> + +h<p 2 ) + 9 2 (m<pF - + hF<p 2 - hipip<j^ . (20) 

The potential energy density term will be as follows: 

[W + h.c.] F = (nKpF - ^mtptp + hF(f> 2 - hi/;ip<pj + h.c. 



(21) 



Kinetic and potential terms averaged over 6 and 9 and emergence of soft SUSY 
breaking terms: 
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Averaging over the Grassmannian coordinates, we compute L =< C >= I d 2 9d 2 9V(9, 9)C. 
Here C is the usual Super-Lagrangian density: C = &$ + W5^(9) + W^5^(9). Then, using 
Eq.17 and Eq.14 < ^kinetic >, namely the kintetic part of L is found as, 



D 



|£| 2 |0| 2 + £*0F* + £0*F. 



Similarly, the potential energy density averaged over 9 and 9 is given by 
<W + h.c >= J d 2 9d 2 9V(9, 9) (W5 {2) (9) + WU { - 2 \9) 
Using Eq.20 and Eq.14 we find, 



(22) 



(23) 



< W + h.c. > 



C Qm0 2 + ^0 3 ^) + (rri(f)F - ~mif)i/) + hFcj) 2 - h^ipc^ 



h.c. (24) 



Total Lagrangian is then: 



L=< C > = < £ kinc tic > + <W + h.C. > . 



(25) 



Using Eqs.21, 22 and 24 we break L into SUSY invariant and SUSY breaking parts as follows. 



L = < C > 



SUSY 



where, 



and 



< C > 



SUSY 



+[W + h.c\ 



Fi 



The equations of motion of auxiliary fields are then, 

F = -(£*0 + m0* + h(j)* 2 
F* = -(£0* + mcj) + h(j) 2 ) 

Substituting F and F* in L one finds, 



C ( \rncp 2 + ~/^ 3 



h.c. 



(26) 

(27) 
(28) 



and 



(29) 



L = L 



On-shcll-SUSY + -^soft, 



(30) 
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where -^on-shcii-susY is the usual on shell SUSY invariant Lagrangian for interacting Wess- 
Zumino model and is given by, 



^On-shcll-SUSY 



2 A. \ 

mh\4>\ <p + -rmpip + hiptp^J 



m 



2n 2 



h.c. 



and L soft is given by, 



-L 



soft 



-F mc 
2^ v 



2 



+ h.c. 



(31) 



(32) 



We remind that a negative sign in the left hand side of Eq.32 appears simply because of 
considering a positive sign before < W + h.c. > in Eq.25 while writing the total Lagrangian. 
We note that it is only the superpotential term in the original theory that leads to soft- 
breaking terms in the resulting Lagrangian (m — > £*m and h — > 2£*h going from W to I/ so f t ). 
Presence of a vector field will not lead to any soft SUSY breaking term. This can easily be 
seen by considering a vector superfield in Wess-Zumino gauge. 



MSSM: 

In MSSM, as mentioned in Ref. [8] the superpotential term [iH u H,i will lead to ^*ixH u Hd 
and terms like y up QU c H u will lead to 2£*y up QU c H u as soft SUSY breaking terms. Here, 
fields with tildes denotes the scalar component of the corresponding chiral superfields. One 
further obtains a gaugino mass term ^-£;A«A«. Thus, one finds a universal gaugino mass 
m i/2 — a bilinear Higgs soft parameter = £*, a universal trilinear soft parameter 
A = 2£*. With no resulting scalar mass term, one has the universal scalar mass parameter 
m = 0. These are the input quantities to be given at a scale A. Low energy spectra is then 
found via RG evoplutions. Ref. [8] considered £ and A as the input quantities and considered 
M G < A < M P . 
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